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SECTION  1 


SCOPE 


1.1  IDENTIFICATION 

Volume  I,  Submarine  OMEGA  Computer  Program  Performance  Specification,  defines 
the  functional  requirements  for  the  Submarine  CKEGA  Computer  Program  which  is 
used  by  the  AN/ARN-99  OMEGA  Navigation  Set.  The  Navigation  Set  and  the  OMEGA 
program  together  comprise  the  Submarine  OMEGA  Navigation  System.  The  tape 
which  defines  the  computer  program  is  entitled  AN/BRN-7  Navigation  Program, 

Volume  II,  Submarine  OMEGA  Computer  Program  Design  Specification,  allocates 
the  functional  requirements  of  Volume  I to  the  computer  routine  and  sub- 
program level. 


This  volume  describes  the  subprogram  designated  as  common  subroutines,  which 
has  no  abbreviation  in  the  listing. 


1.2  COMMON  SUBROUTINE  TASKS 


a)  Compute  Theta  1 

Central  angle  between  station  and  submarine 


b)  Theta  c 

Phase  prediction,  station  to  submarine 

c)  D PHI  KI 

Base  to  station  correction  term 


d)  Position  matrix  computation 

Calculation  of  R. . matrix 
ij 

e)  R,  . rotate 

ij 

Rotate  a 3x3  matrix  through  two  small  angles 

f)  Resolve 

Resolve  submarine  speed  to  component  velocities 

g)  Bearing 

Calculates  bearing  between  system  axes  and  destination 

h)  Vector  by  matrix  multiply 

A 3x3  matrix  by  a 1x3  vector 


i)  Cross  Product 
Vector  product 
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j)  DOT  Product 
Scaler  product 

k)  Square  Root 

Newton's  method  applied  to  square  root 

l)  Sine  Cosine 

m)  Arctangent 

n)  ArcSine 

o)  ArcCosine 
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SECTION  2 

APPLICABLE  DOCUMENTS 


a)  Submarine  CMEGA  Computer  Program  Performance  Specification  (Volume  I 
of  the  Submarine  OMEGA  Computer  Program  Specification). 

b)  Submarine  OMEGA  Computer  Program  Design  Specification  (Volume  II  of 
the  Submarine  C^EGA  Computer  Program  Specification). 

c)  NORT  68-66,  NAP70  User's  Manual,  July,  1968. 

d)  NORT  68-115A,  Detailed  Description  of  NDC-1070  Computer  Instructions, 
Revision  A,  February,  1970. 

e)  NORT  69-87A,  NDC-1070  Flow  Chart  Program,  User's  Manual. 
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SECTION  3 
REQUIREMENTS 


In  order  to  understand  the  program  description  contained  in  the  following 
pages,  it  is  necessary  that  the  reader  will  have  become  familiar  with  the 
associated  functional  requirements  found  in  Volume  I,  Performance  Specification, 
and  with  the  subprogram  allocation  found  in  Volume  II,  Design  Specification, 

3 .1  DETAILED  DESCRIPTION 


3.1.1  Reference  Labels  to  Flow  Diagrams 

The  code  used  to  reference  the  particular  block  in  the  flow  diagrams. 

Section  3.2,  is  as  follows:  The  first  number,  preceded  by  a p,  is  the 
page  number  found  in  the  upper  right  corner  of  the  diagrams.  This  will  be 
followed  by  a slash  sign  (/)  to  separate  the  page  number  from  the  block 
designator.  The  designator  will  either  be  a mnemonic  label  (e.g,,  TEST 
SYNC),  a local  label  indicated  by  a dollar  sign  ($) , or  an  integer.  The 
two  types  of  labels  reference  the  particular  information  block,  on  the 
given  page,  to  which  the  label  is  attached.  The  integer  number,  n,  means 
that  the  referenced  block  is  the  n^^  block  from  the  top  of  the  page;  p8/3 
would  refer  to  page  8 and  the  third  information  designation  down. 

The  label  pl/$2+3  refers  to  page  1,  and  the  3rd  information  block  after 
the  label  $ 2.  p2/7,8,9  refers  to  page  2 and  the  7th,  8th  and  9th  blocks. 

3.1.2  Flow  Diagram  Description 

3)  Compute  Theta  1 (Page  1);  Computes  the  earth  central  angle  between  a 

specified  station  and  the  submarine.  The  arguments  consist  of  the  station 
number,  with  which  the  station  coordinates  can  be  found  by  table  look-up, 
and  a pointer  to  indicate  memory  address  of  the  submarine  position  vector 
to  be  used. 

pl/Compute  Theta, I: 

The  R1  vector  is  converted  to  geocentric  coordinates  by  multiplying  r 

2 2 
by  a /b 

where  a = earth  polar  radius 

b = earth  equatorial  radius 
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£1^2: 

The  Cross  product  and  DOT  product  subroutines  are  used  to  compute 
Sine  = y(R^  X . (Rj^  X S) 


Cos0  * R1  • S 


9 = tAN*^  ^ 

C0S9 

h)  Theta  C (Page  2):  This  routine  computes  the  predicted  phase  from  a station 
to  the  submarine  for  any  frequency.  The  subroutine  D PHI  KI  is  used  to 
compensate  for  the  difference  in  time  between  reception  of  signals  from 
the  station  as  opposed  to  that  of  the  base  station.  The  arguments  are 
0^,  frequency  and  station  number. 

p2 /THETA  C: 

0j^  is  converted  to  cycles  and  rescaled. 


p2/2: 


Obtain  02 > 9^  from  Propagation  Prediction  and  add  to  0^. 


£2/4; 


Obtain  A0  correction  from  D PHI  KI . 


sill- 

9 = 9.  + 0.  + 9-  + A0. 

C L A 5 

c)  D PHI  KI  (Pages  3,  4):  Computes  the  correction  term  that  represents  what 
is  expected  to  happen  to  the  base  station  phase  between  the  base  station 
burst  time  and  the  time  of  a non-base  station  burst  for  any . frequency . 
Arguments  are  station  number  and  frequency. 
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p3/D  PHI  KI  through  p4/$3: 


The  base  station  reference  numbers  and  the  non-base  station  reference 
numbers  are  organized  by  station  and  frequency  to  determine  what  station 
is  sending  on  10.2  kHz  when  the  base  station  is  sending  the  frequency 
of  the  input  argument. 

Reference  Number 


Base  (or  station))  (Frequency)  ^ [(Base  or  Station)  \ 
0-7  / ■ \ 0,1,2  ) \ on  10.2  j 

p4/$3+l  through  p4/$3H-5: 


Modulo  8 


This  loop  decrements  the  station  reference  number  by  unity  each  time 
through.  The  burst  collection  time  is  from  a table  which  lists  the 
span  of  each  burst  in  the  pattern.  The  0.4  seconds  between  the  End 
Burst  and  Start  Burst  is  added.  Exit  when  station  reference  number  is 
equal  to  Base  Station  reference  number. 


p4/$4: 


A 0 = — 0 . 

10  base 


d) 


Position  Matrix  Computation  (Page  5):  This  routine  computes  the 


matrix  from  a latitude  and  longitude.  Used  in  the  Navigation  routine. 


e) 


R^.  Rotate  (Page  6):  Using  a small  angle,  here  T2  and  T3: 


'■12^ 


A0  = 


+ T3 

2 

T3 

-T2 

•T3 

-T3^ 

T2T3 

2 

2 

T2 

T2T3 

2 

-T2 

2 

2 

then 


[r]  . [r]  + [ar]  [r] 


f)  Resolve  Subroutine  (Page  7):  The  Resolve  subroutine  is  used  by  velocity 
processing  (Volume  VIII)  and  in  the  calculation  of  sea  current  (Volume  X) 
to  compute  the  uncorrected  system  velocities  and  V^. 


+ Cos  (9_  + 


h ■ V * *A> 


''CT  =‘"  + *A> 


Where  V 


CT 


0. 
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g)  Bearing  (Page  8):  Computes  the  angle  between  the  axis  and  a vector 
in  the  plane  of  as  indicated  below. 


10759 

'^A  ^ Submarine  Heading  HDG 

= Craft  Velocity  Track  - TK 

D *=  Course  to  Selected  Destination 

4^^  = Sea  Current  Direction 

Op  = System  Heading  N to  CCW 

FIGURE  1 SYSTEM  AZIMUTHAL  RELATIONSHIPS 


1 
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h)  Vector  by  Matrix  Multiply  (Page  9):  Calculates  the  multiplication  of  a 
3x3  matrix  by  a vector  V as  follows: 


(V^.  V^.  V3)  x/  c 


11 

*^12 

^13 

‘21 

^22 

^23 

31 

^^32 

r 

33 

(X3,  X^,  X3) 


X.  = 
1 


i)  Cross  Product  (Page  10):  Calculates  the  vector  product  of  two  vectors: 

X = "U  X “V 

= (U2V3  - U3V2)  i + (U3V^  - U^V3)  T + (U^V,  - U^V^)  k 

j)  DOT  Product  (Paxe  11):  Calculates  the  scaler  product  of  two  vectors: 

X = u^'v 

= -f  U2V2  + U3V3 

k)  Square  Root  (Page  12-13):  This  routine  replaces  the  long  fraction,  F, 
in  registers  0 and  1 by  its  approximate  square  root.  No  other  register 
is  altered.  If  F is  negative,  the  square  root  computation  is  not  carried 
out,  and  a negative  result  is  returned.  The  computation  consists  of  the 
following  steps: 

1)  F is  normalized;  i.e.,  transformed  into  G by  M left  shifts,  such 
that  0.5<G  < 1. 

2)  The  first  (single-precision)  approximation,  Xq,  is  obtained  by: 

Xq  = 0.313567  + (0.890194  - 0.204406  * G)  * G 

3)  An  improved  estimate,  X^^ , is  obtained  by  a single-precision 
evaluation  of  Newton's  formula:  X^  = (X^  + G/Xq)/2 

4)  Newton's  formula  is  evaluated  in  double-precision  to  produce  the 
final  estimate,  X^ . 
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5)  The  result  is  obtained  from  X2  by  adjusting  for  the  scaling  in 
step  1),  thus; 

= X2  * 2'^^^ 

-31 

A maximum  absolute  error  less  than  2 is  expected , with  a negative 
bias  caused  by  truncation  occurring  in  the  right-shift  instruction. 

The  execution  time  is  given  by: 


296  + 4N  microseconds  if  N is  odd  (1  through  31). 

318  + 4N  microseconds  if  N is  even  (2  through  30). 
where  N is  the  number  of  leading  binary  zeros  in  F. 

For  F = 0 the  time  is  54  microseconds. 

The  storage  requirements  are  37  words,  plus  access  to  an  ordered  table 
of  powers  of  two, 

1)  Sine  and  Cosine  (Pages  14-16);  These  routines  take  an  angle  (scaled  tt) 
in  registers  0 and  1 and  returns  results  (scaled  2)  as  follows: 


Entry 


Result 


Other  Registers 


SIN 

COS 

SIN-COS 


sine  in  0 and  1 

cosine  in  0 and  1 

sine  in  2 and  3 and 
cosine  in  0 and  1 


Unchanged 
Unchanged 
Reassigned  +2 


The  sine  computation  takes  the  following  steps; 

1)  The  angle  X is  rescaled  'n^/2  by  the  transformation 

Y * 2X  if  - 1/2  < X <1/2 

Y = 2 (1-X)  if  1/2<X<1 

Y = 2 (-1-X)  if  -1<  X<-l/2 

2)  If  I Y I < 62/109,  the  sine  is  obtained  from  the  formula; 

Sin  X * Y * (Y^) 

(where  S^(t)  is  a polynomial  in  t of  degree  3). 

3)  If  I Y |>  62/109,  the  sine  is  obtained  by  computing 

|sin  X I - C^(Z^) 
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and  adjusting  for  the  sign  of  X (where  Z = 1 - | Y | , and 
C^(t)  is  a polynomial  in  t of  degree  3). 

The  cosine  is  obtained  by  computing  sin  7i(X  + 1/2). 

“8 

The  maximum  relative  error  is  approximately  10  , 

The  maximum  execution  times  for  the  different  entries 


or  2 
are : 


•27. 


SIN  204  microseconds  ] 

COS  208  microseconds  { 

SIN-COS  446  microseconds  ! 

) 

The  storage  requirement  is  58  words.  j 

The  coefficients  for  the  polynomials  S^  and  C^  were  obtained  from  "Computer  | 

Approximations"  (Hart  et.  al.,  1967)  by  making  adjustments  for  scaling  to 

those  given  under  sections  SIN  3180  and  COS  3700.  j 

S^(t)  = 0.785398156180  - 0, 322981354578t 

+ 0.039835751A36t^  - 0.002289069985t^ 

CjCt)  = 0.499999995214  - 0. 1542 12 3732 34t 

+ 0.0079258969611^  - 0.0001603038041^ 


m)  Arctangent  (Pages  17,  18);  This  routine  computes  the  angle  A from  two 
(2)  arguments  S and  C of  the  form  K sin  A (in  registers  2 and  3)  and  K 
cos  A (in  registers  0 and  1).  The  result,  scaled  'ir,  is  returned  in 
registers  0 and  1,  and  two  (2)  registers  are  pruned. 

Use  is  made  of  the  identify 


arctan 


sin  t 
cos  t 


) 


TT 

4 


+ 


arctan 


sin  t - cos  t \ 
sin  t + cos  t / 


The  computation  steps  are 


1)  Transform  S and  C to  a single  argument,  where 


X 


I al  - Id 

I ^ + Tcl 


2) 


Compule  A = Q/8  + X * (X^)  * sign  (S)  * 
where  T^(t)  is  a polynomial  in  t of  degree 
bits'*  Q are  derived  as  follows: 


sign  (C) 

7,  and  the 


"quadrant 
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_9. 

001 

oil 

101 

111 


s 

+ 

+ 


+ 


-7 


or  2 


-23 


The  maximum  relative  error  is  less  than  10 
The  maximum  execution  time  is  518  microseconds. 
The  storage  requirement  is  51  words. 


The  coefficients  of  the  polynomial  T^  were  obtained  from  "Computer 
Approximations"  (Hart  et.  al.,  1967)  by  dividing  those  given  under 
section  ARCTN  4991  by  tt,  and  are  as  follows: 


T^(t)  = 0.318309854667  - 0. 10609902 1806t 


+ 0.063565605401t 


+ 0.031558174030t 


0.044625393671t“ 


0.018935106721t- 


+ 0.007719779427t^  - 0.001493916081t^ 

n)  Arcsine  (Page  19):  This  routine  computes  arcsine  by  calculating  the 
cosine  then  using  the  arctangent  subroutine. 

o)  Arccos  (Page  20):  Computes  the  arccosine  by  calculating  the  sine  then 
using  the  arctangent  subroutine. 


3.2  FLOW  DIAGRAMS 

The  common  subroutine  flow  diagrams  are  presented  on  the  following  pages. 
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QOOl 


• C(m;TC  TICTAI 

• 

• this  ROUTIHC  COMPUTCS  T»C  CCNTftAL  ANOLC  KTftCCN  ANY  STATION  MO  TIC 

• CRAFT.  THt  ARGUMENTS  CONSIST  OF  T»C  STATION  NUMKR  A»0  A ROINTCR  TO 

• THE  POSITION  VECTOR  OF  THE  CRAFT. 


1 


J 

^AQC  0002 


* T>CTA  C 

• 

• this  pout  INC  MILL  COMPUTC  THE  PCOICTEO  PHASE  TROH  A STATION  TO  TIC 

• CPArT  TOP  ANY  PPEQUCNCY . IN  PHASE  DirrCPCNCE  NAVIGATION  A COPPECTION 

• TERM  IS  ACX)C0.  THIS  TERH  IS  T«  TIME  BETICEN  BASE  STATION  BURSTS  AM) 

• THIS  STATION'S  BURSTS  TIMES  THE  PHASE  ESTIMATE  IN  THE  BASE  STATION 

• TRACKING  riLTER  FOP  THIS  FPEOUENCY.  THE  ARGUMENTS  CONSIST  OF  TICTAI» 

• FREQUENCY  AND  STATION  NUMBER. 


BBK 


n 
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0 PHI  K| 


this  poutinc  computes  a coprcction  tcpm  that  pcppcscnts  iaut  is 

exPCCTEO  TO  HAPPEN  TO  THE  OASC  STATION  PHASE  eCTUCCN  THE  BASE  STATION 
BURST  TIME  AND  the  TIME  or  A NON  BASE  STATION  BURST  FOR  ANY  FREQUENCY. 
THIS  ROUTINE  COMPUTES  THE  TIME  DIFFERENCE  AND  USES  THE  PHASE  ESTlHATC 
IN  THE  BASE  STATION  TRACKING  FILTER  AS  T>C  EXPECTED  CHANGE  IN  PHASE  OT 
THE  BASE  STATION  DURING  A TEN  SECOND  PERIOD.  THE  ARGUMENTS  CONSIST  OT 
THE  station  number  AND  THE  FREQUENCY. 


D PHI  Kt 


SET  DELTA  T • 0 


USE  THE  BASE  STATION  NUI«ER  AM) 
FREQ  TO  DETERMINE  UHAT  STATION 
IS  SENDING  10.2  WHEN  T«  BASE  IS 
SENDING  THIS  FREQ.  (BASE  REFI 


□ S BASE  REF  IN  RANGE  OF  0 TO  7T 


PUT  BASE  REF  IN  RANGE  OF  0 TO  7 


USE  THE  STATION  NUMBER  AND  FREQ 
TO  DETERMINE  HHAT  STATION  IS 
SENDING  10.2  HHEN  THIS  STATION 
IS  SENDING  THIS  FREQ  ISTA  REFI 


S STA  REF  IN  RANGE  OF  0 TO  77 


PUT  STA  REF  IN  RANGE  OF  0 TO  7 


• IF  TIC  STATION  NUMBER  IS  TIC 

• BASE  STATION  NUMBER  THE  TEST 

• BELOM  HILL  CAUSE  THE  CORR- 

• ECTION  TO  BE  ZERO 
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PACe  0009 


• POSITION  MATRIX  COMMTATION 

• this  routines  computes  a nitc  CLeHENT  position  matrix  with  t»c  NuracR 

• 3 AXIS  POINTING  NORTH.  THE  ARCUHENTS  CONSIST  OF  THE  LATITUDE  ATO 

• LONGITUDE 


J 


PAOC  0006 


• 6U  OOTATC 

• 

• THIS  PROGRAM  WILL  ROTATC  A 3X3  MATRIX  (Rl 

• THROUGH  ThC  i ANGLES  SPCCiriCO 


A ROTATE  RIJS 

OrvEN  THE  a ANGLES  TZ  ATO  ft 
THEN 


SO 

T3 

-T2 

theta  • -T3 

-T3*T2/2 

T2*T3/2 

T2 

T2*T3/2 

i*CRE  SO  - fT2*Ta  ♦ T3«T3I  / 2 
A»0 

RfNEH)  • RCOLOl  ♦ RCOLO)*T«TA 
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PAGC  0007 


• RCSOLVe  SUBROUTIIC 

• 

• the  resolve  subroutine  is  used  by 

• velocity  PROCCSSINO  to 

• COMPUTE  va  AND  V3  FROM  SHIPS  SPEED  MlTH 

• CROSS  TRACK  VELOCITY  EQUAL  TO  ZERO. 


i RESOLVE 

RESOLVE  ALONG  TRACK  VEL  (REG  21 
AND  CROSS  TRACK  VEL  (REG  31 
THROUGH  ANGLE  IN  REG  O.l^-THETA  P 
TO  FORM  V2.V3  (RESCALED  BY  Bll 
LEAVE  RESULT  IN  REG  2,3  AFTER 
PRUNING  2 REGISTERS 
V2  • VAT*S1N  ♦ VCT*COS 

V3  ■ VAT*CQS  - VCT»SIN 

B8K 


j 

i! 
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>>Acc  oooe 


• BCARINS 

• 

• THIS  ROUTINE  MILL  COMPUTE  THE  ANGLE  BETWEEN  THE  R3  AXIS  AN>  A VECTON 

• IN  THE  PLANE  OT  Re.RS  POINTING  AT  A FIXED  POSITION.  THE  ARGUMENTS 

• CONSIST  or  A POINTER  TO  THE  R^  VECTOR  OT  THE  RIJ  MATRIX  AND  A POINTER 

• TO  A VECTOR  THAT  OEFINES  THE  FIXED  POSITION. 


J 


PAOC  0009 


• VCCTOR  BY  MATRIX  MULTIPLY 

• 

• THIS  ROUTINE  multiplies  A 3X3  MATRIX  (R)  BV  A VECTOR  U TO  GET  A NEW 

• VECTOR  X.  the  arguments  consist  or  POINTERS  TO  R.  U ATO  X. 


« 


r' 


J 

• DOT  phoouct 

• this  ROUTIIC  rORHS  T»c  DOT  PROCXJCT  QT  THE  THO  VECTORS  U AND  V 

1 REGISTER  DOT 

U . V - SUH  Ufn*Vltl.  t«t  TO  3 

aSK 


PMX  OOlS 
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C0MPI,1MCNT  RESULT  2(1-  X)  • Y 
ir  X IN  2N0  OR  -8(1  ♦ XI  • V 
ir  X IN  3RO  QUADRANT 


PAOC  0019 
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PAOC  0017 


• ARCTANGCNT 

• this  routine  mill  compute  t»c  arctangent  for  a sine  cosine  pair, 

• T«  SINE  AND  COSINE  ARCUHENTS  MUST  HAVE  THE  SAME  SCALE  FACTOR  A»«>  T»C 

• RESU.T1NG  ANCLE  IS  SCALED  PI. 

• 


i ATAN  ENTRY 

ADJUST  RIS  STACK 


« ENTER  HERE  HHEN  NOT  USINO  PIN  CNTRV 
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coMPure  cos 

• <1  - SIN»«2)»«.5 

SCT  cos  - 0 IF  SIN  none  than  i 

COMPUTE  angle 
- ADCTAN  SIN/COS 

(EXIT  ASIN  FROM  AffCTANt 


ROP 


Electronics  Division 


NORTH 


3.3  CC»1PUTER  subprogram  ENVIRONMENT 


3.3.1  Common  Subroutine  - Tables 

a ) D PHI  KI  Routine: 

Station  Burst  Time  Table:  This  table  is  used  to  determine  the  burst 
time  for  each  station.  The  table  is  defined  in  detail  in  the  flow 
charts  (page  14  of  CMEGA  Processing)  and  in  the  listing.  The  resolution 
for  the  burst  time  is  0.005  second.  The  official  name  of  the  table  is 
STATION  BURST  TIMES. 

b)  TOETA  C Routine; 

Lambda  Table;  This  table  is  used  to  convert  radians  to  cycles  at  the 
appropriate  frequency.  It  is  defined  in  the  listing.  The  official  name 
of  this  table  is  LAMBDA  TABLE. 

c ) THETAl  Routine; 

Station  Locations  Table;  This  table  contains  the  locations  of  all 
existing  CMEGA  transmitting  stations.  Each  location  is  specified  as 
a three  element  geocentric  position  vector.  The  first  entry  is  for 
station  A.  The  table  is  defined  in  detail  in  the  listing.  The  official 
name  of  this  table  is  STATION  VECTOR  TABLE. 


3.3.2  Temporary  Storage 

All  temporary  storage  is  in  the  R15  push-down  stack. 


3.3.3  Input/Output  Formats 
Not  applicable. 
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3.3.4  Required  System  Library  Subroutines 

All  subroutines  listed  below  are  found  in  this  volume. 


SUBROUTINES 

CALLING  LABEL 

FLOW  DIAGRAM 

ARCCOSINE 

ARCCOS 

p20 

ARCSINE 

ARCSIN 

Pl9 

ARCTANGENT 

ATAN 

pl7  pl9/l 

pl/3  p20/2 

P8/1 

BEARING 

BEARING 

p8 

COMPUTE  THETA 1 

CCMPUTE  THETA  1 

pi 

CROSS  PRODUCT 

REGISTER  CROSS 

DOT  PRODUCT 

REGISTER  DOT 

pl/2  p8/l 

pll 

D PHI  KI 

D PHI  KI 

POSITION  MATRIX 
CCMPUTATION 

PANEL  MAIN  $8 

P5 

RESOLVE 

RESOLVE 

P7 

RIJ  ROTATE 

ROTATE  RIJ'S 

p6 

SIN-COS 

SINCOS 

SQUARE  ROOT 

SQUARE  ROOT 

Pl2 

THETAC 

THETAC 

p2 

VECTOR  BY  MATRIX 
MULTIPLY 

MATRIX  30 

P9 
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